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Abstract. We presented here the orbital parameters for five visual binary 
stars calculated by using the new method which we named Sector Grid Search. 
Orbital parameters were obtained for the following stars: WDS 00152+2722 = 
ADS 195, WDS 02202+2949 = ADS 1780, WDS 11550-5606 = HIP 58106. 
WDS 16256-2327 = ADS 10049 and WDS 16256-2327 = ADS 10045. In 
addition, their masses, dynamical parallaxes and ephemerides were calculated 
as well. 
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1. INTRODUCTION 

The study of binary stars is a very old branch of astronomy which still plays 
an important role in a present astronomy. It has a number of aims. Binary stars 
have been studied for decades to measure accurate the stellar masses, to test the 
evolutionary models and the star formation theories. 

Although the Washington Double Star Catalog (WDS) (Mason et al. 2006) 
contains a large number of double and multiple systems (over 84000), according to 
the Sixth Catalog of Orbits of Visual Binary Stars (Hartkopf & Mason 2006) only 
for about 1900 systems the orbital parameters have been calculated. It is because 
available data for many systems are not sufficient for an accurate determination 
of the orbital parameters. It is caused by two main reasons. The number of 
observers is too small compared with the number of systems and/or the orbital 
period in some systems is too long compared with a time span since the beginning 
of systematic measurements. As a result researchers have been pushing to estimate 
system parameters from ever shorter arcs. The literature on the determination of 
system parameters of visual binary is a very numerous. Despite such a large 
number of methods, no one method proposed so far is adequate in the case of 
a short arc. We presented here a fully automated method for determination of 
systems parameters from a short arc, which we named Sector Grig Search (SGS). 
Furthermore, we also presented the systems parameters of five visual binary stars 
calculated by using SGS method. The orbits of the stars WDS 00152+2722, WDS 
02202+2949 and WDS 16256-2327 (pair AB) were previously announced by IAU 
Commission 26 (2007, Inf. Circ. 162). 

The dynamical parallaxes and the individual masses were calculated for stars 
which belong to the Main Sequence according to the method proposed by Angelov 
(1993). Trigonometric parallaxes published in the Hipparcos and Tycho Cata- 
logues (ESA 1997) were used in order to calculate the total masses of the systems. 
The values of the stellar masses as a function of the spectral type, used as refer- 
ence, have been taken from Schmidt-Kaler (1982). Visual magnitudes and spectral 
types, presented in Table 3. were taken from WDS catalogue. 
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2. METHOD 

Among the methods presented so far one can find several efficient methods. 
We will mention here only a few of them which are widely used. 

The best known method for determination of system parameters is that of 
Thicle-Innes-van den Bos (Thiele 1883; van den Bos 1926, 1932; Dommanget 
1981). Yet it cannot be expected to handle all possible cases. It requires the 
knowledge of three normal places and an area constant. Docobo (1985) proposed 
an analytical method which does not require the knowledge of the area constant. 
The method is based on a mapping from the interval (0,2n) into the family of 
Keplerian orbits whose apparent orbits pass through three base points. These 
points are taken either coinciding with the most reliable measurements or belong 
to the areas with a maximum of observational evidence in their favor, but do not 
necessary coincide with actual measurements. The orbits which fits the best to 
all known measurements is chosen from all generated possible ones. In addition, 
it has an option allowing the use the radial velocities, being thus useful for the 
inclinations close to 90° and in the case of a short arc as well. 

The numerical methods, based on minimization procedure, which take advan- 
tage from computer's power were proposed by Eichhorn and Xu (1990) and Pour- 
baix (1994). The method proposed by Hartkopf et al.(1989) also takes advantage 
from computer power but with a different approach. The orbital parameters are 
calculated according to the "three-dimensional adaptive grid search" and it re- 
quires preliminary knowledge of the period P, epoch of periastron passage T and 
eccentricity e. 

The most recently Olevic & Cvetkovic (2004) proposed a method suitable for 
application in the case of short arc which is based on using supplementary, fictive, 
measurements. In our previous papers (e.g. Novakovic 2007) we used this method 
in order to calculate the orbital parameters and we obtained acceptable results 
when deal with a short arc. Although we obtained good results with this method, 
they depend on the choice of the fictive measurements. 

We presented here a simple and fully automated method which can give a 
reliable orbit solution even in the case of a short arc and does not need any 
subjective initial parameters such as fictive measurements. It consists of several 
main steps. 

The first step in our procedure is to assign the appropriate weights to all mea- 
surements according to the weighting rules described in Hartkopf et al. (1989,2001) 
and rejecting all measurements with errors which do not satisfied 3cr condition. 
The next step is to divide observed arc in five equal intervals and determining min- 
imum and maximum value of angle separation for every interval (p m in, Pmax)i=i,5 
as well as an average value of position angle for every interval 0j=i,5. El Given the 
set of the five position angles 0j=i,5 and the five intervals of the angle separations 
{pmin, Pmax) i=i,5 we can proceed further. As we know a general equation of the 
conic section is of the form 



1 ln order to include the possible measurements errors it is necessary to correct the values of 
Pmin and pmax (correction makes interval a little bit wider). This correction should takes into 
account a maximum allowed error of 3o\ 
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(1) 



so the curve has been determined if we have determined the five coefficients 
O20:ao2,aii,ai0iaoi- If the positions of the five points were given the coefficients 
can be determined by solving system of the five linear equations. We apply an 
iterative procedure beginning with a p m in and increasing the angle separation p 
for step size as many times as we need to reach p m ax- For every combination 
of the five position angles and the five angle separations the five coefficients arc 
determined by using LU decomposition described in Press et al. (1992). From 
these five coefficients we can calculate the five orbital parameters (a, e, i, f2, oj) and 
two dynamical parameters (P, T) can be calculated from the Keplerian equation 



In this formula, Wj denotes the weight of the j th measurement, x ,y denotes 
the observed positions and x c ,y c denotes the calculated positions of the companion. 
The minimization is done by applying Powell's method (Press et al. 1992), but 
with just three degrees of freedom (P, T, a) in order to decrease probability to reach 
a local minimum which number growing as exp(iV) (V-the number of degrees of 
freedom) . For each step in the grid search a set of orbital parameters is calculated 
and a solution with minimum residuals is accepted as a final solution. 



In order to test efficiency of SGS method we determined the orbital parameters 
of two visual binary stars which have had orbital parameters already calculated 
and very reliable. For this test we selected the star WDS 02022+3643 (HD 12376) 
from the Catalog of Orbits and Ephemerides of Visual Double Stars (Docobo et 
al. 2001) and the star 16413+3136 (HD 150680) from the Sixth Catalog of Orbits 
of Visual Binary Stars. In the case of the star WDS 02022+3643 test was made 
using speckle measurements only, but in the case of the star WDS 16413+3136 
test results were obtained using mixed data (visual, speckle, photographic, CCD). 

The star WDS 02022+3643 (pair AB) has a very good orbit calculated by 
Hartkopf et al. (2000). As we wanted to show that method presented here is 
efficient in the case of a short arc we did not use all measurements available for this 
star but only a set of 7 measurements from 1994.7086 to 1999.8856 which covered 
an arc of 59°. Similarly, we calculated orbital parameters for binary star WDS 
16413+3136 using only a part of available measurements. For the purpose of this 
study we chose 38 measurements from 1989.3931 to 1999.61 which covered an arc 
of « 67°. 

Test results are presented in a Table 1. and they show that by using SGS 
method good results can be obtained even in the case of a short arc. 



— (t j -T)=M(6 j ,a,e,i,n,oj). 
The last step is minimization of function D defined as 



(2) 




(3) 



2.1. Test 
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Table 1: Test results 



Element 


SGS method 


Hartkopf ct al. 2000 


SGS method 


Sodcrhjelm 1999 


P[yr] 


12.90 


12.94 


36.76 


34.45 


T 


1988.82 


1989.06 


1966.1 


1967.7 


«["] 


0.155 


0.150 


1.40 


1.33 


e 


0.508 


0.404 


0.40 


0.46 




68.9 


67.0 


127.4 


131.0 




195.2 


191.4 


53.7 


50.0 




289.0 


295.1 


107.4 


111. 



3. RESULTS AND COMMENTS 



In this section we presented results and comments on individual objects. 

Table 2. presents the corresponding numerical values for the orbital parameters 
(epoch J2000) and their estimated formal errors as well as the identifications of 
the stars in several widely-used catalogs. 

Table 3. gives the astrophysical quantities for both components: visual mag- 
nitudes, spectral types, absolute visual magnitudes, masses and, in the last two 
columns, the calculated dynamical parallax and the Hipparcos trigonometric par- 
allax. 

Table 4. gives predicted ephemerides for these systems for the period 2008- 
2012. 

The orbits are illustrated in Figures 1-5. The solid curves represent the newly 
determined orbital parameters, while the dashed curves represent previously pub- 
lished orbital parameters. The solid lines indicate the line of nodes. All measure- 
ments are connected with their predicted positions on the new orbit by "O-C" 
lines. The interferometric measurements are represented of each plot by fulled 
circles "•" and all others measurements (visual, photographic) are represented by 
plus " +" . The direction of motion is indicated on the north-east orientation in 
the lower right of each plot and "+" indicate position of primary star. 

WDS 00152+2722 = J 868. It was discovered in 1907.78 by Barton (1926) 
and 15 measurements have been performed till a year 2002. Recently, an additional 
measurement was carried out by our double star group (Cvetkovic et al. 2007). 
Our orbital solution is the first one for this binary star and it was calculated by 
using set of all 16 measurements. As data about this star are very poor we do not 
know very much about it. Under the assumption that its components belong to 
the Main Sequence we estimated the total mass of system Mtot ~ 1.0 M®. 

WDS 02202+2949 = A 961. This binary star was discovered in 1905.87 
by Aitken (1906) at Lick Observatory. Previous orbit of this binary have been 
calculated by Heintz (1969). In order to demonstrate improvement of previous 
orbit obtained by Heintz we calculated root mean square (rms) of residual values^. 
The obtained rms values are Ad = 2.9°, Ap = 0'.'039 (this work), and AO = 
2.9°, Ap — 0'.'042 (Heintz's orbit) what indicates that our orbit solution does 
improve separation but gives statistically same result for position angle. Our 
orbital parameters and the Hipparcos parallax yield a total mass of system Mtot 



2 When deriving these values corresponding weights have been taken into account 
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ss 3.0 M.q. This value is slightly larger then expected value for this spectral type 
(on condition that both components belong to the Main Sequence). 

WDS 11550-5606 = HLD 114. This is a wide binary star discovered in 
1882 by Holden (1884). Till the present time 47 measurements have been per- 
formed and they cover an arc of «60°. Our orbital parameters are the first ones 
calculated for this binary and they should be classified as preliminary. The new 
measurements are necessary to obtain more reliable orbit solution and to under- 
stand significant difference between dynamical parallax (this work) and Hipparcos 
trigonometric parallax. In general, this disagreement could be due to the uncorrect 
orbital solution and/or because at least one of the components does not belong to 
the Main Sequence. The spectral type of primary star (G3IV/V) suggested that 
it could lies outside of the Main Sequence. 

WDS 16256-2327 = H 2 19 AB. It was discovered as early as in 1780.34 
by W. Herschel. Though the time interval after the discovery has already exceeded 
200 years, available measurements cover a short arc of » 28°. This pair is included 
in the Catalog of Rectilinear Elements (Hartkopf et al. 2006) indicated that it is 
not a physical. Despite this indication, we calculated system parameters because 
we think that is not clear yet whether this pair is physical or optical. The rms 
values are AO = 1.9°, Ap = 07147 (this work), and AO = 1.8°, Ap = 07148 
(linear fit) what is statistically indistinguishable results. Our orbital parameters 
and Hipparcos parallax yield a total mass of the system M. tot ~ 23.6 A4q. As one 
can expect near IO.OtVJ© (Schmidt-Kaler 1982) for stars of B2V spectral type and 
slightly over IQ.QMq for B2IV spectral type, the obtained value of total mass of 
the system is reasonable. 

WDS 16256-2327 = BU 1115 DE. This binary star belong to the same 
multiple system as the previous one, but it was discovered by Burnham (1894) in 
1889.39 at Lick Observatory. The orbital parameters presented here are the first 
for this binary. The dynamical parallax is n — 0700769, in very good agreement 
with the Hipparcos parallax tthip = 070 0733, but the total mass of the system 
{M-tot ~ 5.8 Mq with Hipparcos parallax), calculated with our orbital parameters 
and either Hipparcos or dynamical parallax, is significantly smaller then expected 
for B5V spectral type. As we know, combination of the uncertainty of the semi- 
major axis, of the orbital period and that of the parallax is easily propagated in 
the standard deviation of the total mass as, 




According to equation (j4|) the error of total mass is au = 4.0 M.q- Taking into 
account this error, we obtained 9.8 Mq as the upper limit for the total mass of 
this system. 
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Table 2. Orbital parameters (J2000). 



Name 
WDS 



ADS 
BDS 



HD 
HIP 



P[yr] 



J 868 
00152+2722 



195 



1089. 1725. 
±153. ±147. 



8.48 
±1.06 



0.367 75.2 
±0.370 ±2.5 



68.9 352.8 
±3.1 ±19.4 



A 961 1780 14394 143.1 2038.1 0.294 0.571 145.4 150.0 15.6 

02202±2949 12885 10892 ±5.6 ±5.7 ±0.057 ±0.095 ±3.9 ±5.8 ±8.5 



HLD 114 
11550-5606 



103493 930. 2039. 8.02 0.707 97.7 167.4 59.7 
58106 ±98. ±91. ±0.90 ±0.109 ±1.5 ±1.3 ±7.4 



H 2 19 AB 10049 147933 2398. 2327. 4.25 0.675 135.3 77.5 226.1 

16256-2327 7613 80473 ±326. ±343. ±0.79 ±0.322 ±6.9 ±13.5 ±15.3 

BU 1115 DE 10045 147888 675.5 2008.6 1.01 0.707 134.8 152.7 260.4 

16256-2327 7609 80461 ±32.5 ±34.2 ±0.15 ±0.112 ±2.7 ±4.6 ±1.1 



Table 3. Dynamical elements. 



WDS 



niA-niB 



Sp. 



Ma M b Ma Mb Tr dyn tvhip ± a nHI 



00152+2722 12.42 - 12.62 - 11.91 12.21 0.54 0.50 79.06 

02202+2949 9.31 - 8.94 F5 3.46 3.09 1.33 1.41 7.69 7.48 ± 2.59 

11550-5606 7.36 - 7.81 G3IV/V 6.18 6.63 0.88 0.82 70.50 32.33 ± 1.44 

16256-2327(AB) 5.07 5.74 B2IV B2V ----- 8.27 ± 1.18 

16256-2327(DE) 7.04 8.65 B5V 1.47 3.08 3.06 1.97 7.69 7.33 ± 1.37 



Table 4. Ephcmerides. 

2008 2009 2010 2011 2012 

wds Desig. e P e P e p e p e P 

m(2ooo) n n n n n n n n n n 

00152+2722 228.3 5.817 228.5 5.855 228.6 5.892 228.8 5.929 229.0 5.967 

02202+2949 277.4 0.301 275.7 0.294 274.0 0.287 272.1 0.280 274.2 0.273 

11550-5606 169.7 3.291 169.4 3.264 169.2 3.235 169.0 3.205 168.7 3.172 

16256-2327(AB).. 338.6 2.913 338.5 2.908 338.3 2.903 338.1 2.898 337.9 2.893 

16256-2327(DE).. 259.9 0.214 253.8 0.211 247.7 0.210 241.5 0.210 235.4 0.213 
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Figure 5: BU 1115 DE. 
4. CONCLUSIONS 

We presented here a new method for determination of the orbital parameters 
of visual binary stars. This method yields to the reliable orbit solution even 
in the case of a short arc. On the other hand SGS method is fully automated 
and very easy to use what is specially important when deal with a large amount 
of data. Thus the method will be useful in future astrometric missions such as 
SlM3 and GALA0. In order to improve the orbital parameters of visual binaries 
presented here, new measurements are necessary. The observers should pay the 
more attention to the star BU 1115 DE which will pass through the periastron at 
the middle of 2008. 
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